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Abstract.  Observations  of  polar  stratospheric  clouds  (PSCs)  by  lidar  show  that 
the  clouds  often  contain  solid  particles,  which  are  most  likely  composed  of  nitric 
acid  hydrates.  However,  laboratory  experiments  indicate  that  such  hydrate  particles 
are  not  easily  formed  under  Arctic  synoptic  scale  conditions,  suggesting  that  solid 
PSC  particles  should  be  rather  rare.  Here  we  show  results  from  a  model  study 
indicating  that  mountain-induced  mesoscale  temperature  perturbations  may  be  an 
important  source  of  nitric  acid  hydrate  particles  in  the  Arctic.  Multiple  Arctic 
vortex  trajectories  were  combined  with  a  global  mountain  wave  forecast  model  to 
calculate  the  potential  for  solid  particle  formation  during  December  and  January 
1994/1995.  The  mountain  wave  model  was  used  to  calculate  adiabatic  cooling  over 
several  thousand  ridge  elements.  Nitric  acid  hydrate  particles  were  assumed  to  form 
in  the  mountain  waves  according  to  several  microphysical  mechanisms,  and  were  then 
advected  using  polar  vortex-filling  synoptic  trajectories  to  generate  maps  of  solid 
particle  occurrence.  The  calculations  show  that  mountain  waves  may  be  a  significant 
source  of  PSCs  containing  solid  particles  that  are  observed  on  the  synoptic  scale.  In 
particular,  the  east  coast  of  Greenland,  the  Norwegian  mountains,  and  the  Urals  are 
found  to  be  solid  particle  sources,  with  the  PSCs  often  predicted  to  survive  several 
thousand  kilometers  downstream. 


1.  Introduction 

Polar  stratospheric  clouds  (PSCs)  are  frequently  ob¬ 
served  at  temperatures  above  the  ice  frost  point  in  the 
Arctic  and  are  composed  largely  of  condensed  HNO3 
in  various  forms.  Formation  of  liquid  type  1  PSCs 
due  to  HNO3  uptake  by  sulfuric  acid-water  aerosols 
to  form  large  HNO3-H2SO4-H2O  droplets  is  now  well 
established  [Carslaw  et  al,,  1994;  Drdla  et  a/.,  1994; 
Tabazadeh  et  al.,  1994;  Carslaw  et  a/.,  1997a;  Del  Ne¬ 
gro  et  al.,  1997].  However,  the  mechanisms  by  which 
solid  nitric  acid-containing  particles  form  are  still  un¬ 
clear  [e.g. ,  Peter ,  1997]. 

Lidar  observations  provide  a  means  of  distinguishing 
PSCs  containing  mostly  liquid  particles  from  those  con- 
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taining  solid  particles  [Browell  et  al.,  1990;  Toon  et  al. , 
1990].  Although  there  exists  no  comprehensive  clima¬ 
tology  of  solid  and  liquid  PSC  observations  by  lidar,  the 
many  data  sets  that  have  been  reported  in  the  litera¬ 
ture  indicate  that  solid  particles  are  present  in  a  large 
proportion  of  PSCs  [Browell  et  a/.,  1990;  Beyerle  et  al, 
1994;  Stebel  et  al,  1997;  Nagai  et  al. ,  1997;  Biele  et  al, 
1998;  Muller  et  al,  1998;  Neuber  et  al,  1998]. 

The  phase  of  PSC  particles  has  two  potentially  im¬ 
portant  effects  on  stratospheric  chemistry.  First,  the 
rates  of  heterogeneous  chlorine  activation  reactions  are 
known  to  differ  on  liquids  and  solid  nitric  acid  hy¬ 
drates  [Ravishankara  and  Hanson ,  1996].  However,  re¬ 
cent  modeling  studies  suggest  that  these  differences  in 
reactivity  may  have  a  relatively  minor  effect  on  the  over¬ 
all  rate  of  chlorine  activation  and  subsequent  ozone  de¬ 
pletion  over  the  course  of  an  Arctic  winter  [Carslaw  et 
al,  1997b;  Becker  et  al,  1998].  Of  more  importance 
may  be  the  potential  of  solid  nitric  acid-containing  par¬ 
ticles  to  denitrify  the  polar  winter  stratospheres.  Deni¬ 
trification,  that  is,  the  removal  of  odd-nitrogen  species 
from  the  stratosphere  (principally  HNO3),  promotes  the 
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halogen- catalyzed  destruction  of  ozone.  This  depletion 
in  HNO3  is  a  characteristic  feature  of  the  Antarctic 
during  winter  but  has  also  been  observed  in  the  Arc¬ 
tic  [Arnold  et  al.,  1989;  Hubler  et  al,  1990;  Sugiia  et 
al,  1998].  The  most  likely  cause  of  denitrification  is 
gravitational  settling  of  both  ice  particles  (incorporat¬ 
ing  nitric  acid)  and  large  nitric  acid  hydrate  (NAH)  par¬ 
ticles.  Selective  nucleation  of  a  small  number  of  NAH 
particles  from  an  aerosol  ensemble  allows  these  parti¬ 
cles  to  grow  to  sizes  at  which  they  have  significant  fall 
velocities,  thus  enabling  them  to  remove  nitric  acid  by 
sedimenting  out  of  the  polar  stratosphere.  The  role 
of  selectively  nucleated  NAH  particles  in  denitrification 
remains  to  be  confirmed,  but  it  is  likely  to  be  the  domi¬ 
nant  process  in  the  Arctic  where  temperatures  are  usu¬ 
ally  too  high  to  support  prolonged  ice  existence.  The 
importance  of  solid  nitric  acid-containing  particles  in 
this  respect  makes  their  microphysical  properties  more 
than  of  academic  interest. 

Koop  et  al.  [1997]  have  reviewed  the  processes  leading 
to  phase  transformations  in  PSCs.  Several  mechanisms 
were  identified  that  could  account  for  the  existence  of 
solid  particles  in  Arctic  PSCs.  These  are  summarized 
in  the  following  five  scenarios. 


temperatures  are  low  enough  for  the  initial  formation 
of  ice.  This  is  a  strong  constraint  in  view  of  recent 
work  which  suggests  that  liquid  stratospheric  aerosols 
can  be  supercooled  by  as  much  as  3  K  below  the  ice 
frost  point  [Koop  et  al.,  1998]  or  even  4  K  [Carslaw 
et  al.,  1998a].  This  would  imply  temperatures  as  low 
as  184-185  K  on  the  50  mbar  pressure  level  in  order 
to  initiate  the  ice  formation  required  for  NAH  particle 
nucleation.  Ice-mediated  NAH  particle  formation  may 
be  the  dominant  source  of  solid  type  1  NAH-containing 
PSCs  in  the  cold  Antarctic  stratosphere  where  ice  PSCs 
occur  frequently.  However,  large-scale  ice  occurrence  is 
much  less  common  in  the  Arctic  stratosphere  due  to  the 
higher  temperatures,  so  an  ice  particle  mediated  NAH 
formation  mechanism  operating  on  the  synoptic  scale 
is  unlikely  to  explain  the  large  number  of  observations 
of  PSCs  containing  solid  particles  in  the  Arctic.  In¬ 
deed,  many  Arctic  observations  provide  an  indication 
for  solid  type  1  PSCs  in  air  parcels  that  have  remained 
above  the  ice  frost  point  on  the  synoptic  scale  for  several 
days  [Larsen  et  al. ,  1997].  However,  as  we  show  in  this 
paper,  mesoscale  cooling  events  caused  by  mountain- 
induced  gravity  waves  could  increase  the  likelihood  of 
NAH  formation  due  to  ice  clouds. 


1.1.  Scenario  1 

Nitric  acid  hydrates  could  nucleate  heterogeneously 
on  water  ice  crystals  below  the  ice  frost  point,  lead¬ 
ing  to  the  release  of  NAH  particles  upon  subsequent  ice 
evaporation  (shown  as  scenario  1  in  Figure  1).  Nitric 
acid  hydrate  formation  by  this  mechanism  can  in  prin¬ 
ciple  operate  on  the  synoptic  scale,  provided  synoptic 
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Figure  1.  Particle  phase  transitions  leading  to  the 
formation  of  solid  type  1  PSC  particles  (NAT).  The 
numbers  1-5  refer  to  the  scenarios  described  in  section 
1.  Note  that  the  temperature  axis  is  not  to  scale.  A 
diagram  showing  many  more  particle  phase  transitions, 
from  which  this  figure  was  taken,  is  given  by  Koop  et 
al.  [1997]. 


1.2.  Scenario  2 

Nitric  acid  hydrate  particle  nucleation  might  occur  in 
non-equilibrium  HNO3-H2SO4-H2O  droplets  in  moun¬ 
tain  waves,  without  ice  formation  (see  Figure  1).  This 
mechanism  has  been  shown  to  operate  in  principle, 
based  on  model  simulations  and  laboratory  data  [Tsias 
et  al.,  1997],  but  there  is  as  yet  no  observational  evi¬ 
dence  to  confirm  that  solid  PSC  particles  are  generated 
in  the  stratosphere  by  this  mechanism.  A  reliable  simu¬ 
lation  of  NAH  particle  formation  on  the  global  scale  via 
this  mechanism  is  difficult  since  the  heating  and  cooling 
rates  required  [Tsias  et  al. ,  1997]  are  extreme  and  there¬ 
fore  not  easily  parameterized  for  global  application. 

1.3.  Scenario  3 

Tabazadeh  et  al.  [1995]  and  Tabazadeh  and  Toon 
[1996]  suggested  that  nucleation  of  solid  HNO3/H2O 
phases  might  occur  in  air  parcels  that  have  spent  suf¬ 
ficiently  long  below  the  nitric  acid  trihydrate  (NAT) 
equilibrium  temperature.  They  analyzed  total  aerosol 
volumes  in  a  restricted  aerosol  size  band  as  measured  on 
board  the  ER-2  aircraft,  and  by  assuming  the  particles 
were  in  thermodynamic  equilibrium  showed  that  the  ob¬ 
servations  were  inconsistent  with  NAT.  They  therefore 
postulated  the  existence  of  a  new  phase  with  a  higher 
H20  to  HNO3  ratio  than  either  NAT  or  HNO3-H2SO4- 
H20  droplets,  so  called  type  lc  PSCs.  An  analysis  of 
the  measured  particle  size  distributions  underlying  the 
integrated  particle  volumes  [Peter,  1997]  suggests  an  al¬ 
ternative  explanation  in  terms  of  well  established  solid 
phases.  Nevertheless,  as  already  stated,  the  origin  of 
these  solids  remains  unclear. 
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1.4.  Scenario  4 

Zhang  et  al  [1996]  have  suggested  that  nitric  acid 
hydrates  may  form  on  frozen  sulfuric  acid  aerosols  (sul¬ 
furic  acid  tetrahydrate  (SAT))  after  their  surfaces  have 
become  activated  by  an  initial  deposition  of  NAT  (see 
Figure  1).  The  supersaturation  required  for  NAT  nu- 
cleation  on  preactivated  SAT  particles  corresponds  to 
about  a  2-3  K  cooling  below  the  NAT  equilibrium  tem¬ 
perature.  Such  a  preactivation  mechanism  was  shown 
to  operate  efficiently  in  the  laboratory  and  may  indeed 
operate  in  the  stratosphere. 

1.5.  Scenario  5 

Koop  and  Carslaw  [1996]  have  shown  that  SAT  parti¬ 
cles  can  deliquesce  to  form  HNO3-H2SO4-H2O  droplets 
upon  cooling  to  temperatures  a  few  K  above  the  ice 
frost  point  (typically  3-4  K  below  the  NAT  equilibrium 
temperature,  see  Figure  1).  They  suggested  that  this 
mechanism  might  lead  to  the  heterogeneous  nucleation 
of  NAT.  Recent  laboratory  experiments  by  Iraci  et  al . 
[1998]  appear  to  confirm  the  heterogeneous  nucleation 
of  NAT  on  deliquescing  SAT.  However,  a  similar  study 
[Martin  et  al.,  1998]  indicates  that  the  deliquescence  of 
SAT  may  not  occur  readily  in  the  stratosphere.  This 
NAT  formation  mechanism,  if  effective,  would  com¬ 
pete  with  the  preactivation  mechanism  of  Zhang  et  al. 
[1996].  The  extent  to  which  NAH  particles  are  pro¬ 
duced  by  either  mechanism  4  or  5  depends  on  the  how 
many  of  the  background  aerosols  are  present  as  SAT, 
There  is  presently  little  evidence  for  SAT  particles  in 
the  stratosphere,  which  is  unlikely  to  form  by  homoge¬ 
neous  nucleation  of  liquid  aerosols,  but  rather  via  het¬ 
erogeneous  nucleation  on  either  NAT  or  ice  [Koop  et 
al.,  1997].  Hence,  NAT  formation  via  the  Zhang  et  al. 
[1996]  or  Koop  and  Carslaw  [1996]  mechanisms  may  oc¬ 
cur  in  air  parcels  that  have  already  supported  ice  PSC 
formation.  We  return  to  this  possibility  later. 

All  of  these  mechanisms  have  conceptual  difficulties 
if  they  are  to  account  for  the  presence  of  solid  particles 
in  type  1  PSCs  in  the  Arctic  [Peter,  1997].  The  direct 
formation  of  NAH  PSCs  from  liquid  HNO3-H2SO4-H2O 
aerosols  was  discounted  by  Koop  et  al.  [1997]  on  the  ba¬ 
sis  of  laboratory  experiments.  Scenario  1  requires  tem¬ 
peratures  to  fall  below  the  ice  frost  point,  which  occurs 
rarely  on  the  synoptic  scale;  scenario  2  has  yet  to  be 
confirmed  by  field  observations;  scenario  3  relies  on  the 
difficult  comparison  of  observed  aerosol  volumes  and 
equilibrium  models  to  confirm  a  new  solid  PSC  type; 
scenarios  4  and  5  rely  upon  the  existence  of  SAT  parti¬ 
cles,  which  at  best  form  only  after  ice  PSCs  have  already 
existed. 

In  this  work  we  consider  the  possibility  that  solid 
PSC  particles  are  formed  via  mechanism  1  in  mountain- 
induced  gravity  waves  rather  than  on  the  synoptic  scale. 
The  recent  work  of  Carslaw  et  al.  [1998a]  supports  this 
possibility.  Their  analysis  of  aircraft  lidar  observations 
reveals  that  ice  clouds  that  formed  in  mountain  waves 
over  Norway  acted  as  a  source  of  solid-containing  type  1 


PSCs  that  persisted  downwind  of  the  mountains.  Fur¬ 
thermore,  Carslaw  et  al.  [1998b]  have  calculated  that 
for  the  winter  1994/95,  temperatures  may  have  passed 
below  the  ice  frost  point  in  mountain  waves  as  often  as 
once  every  17  days  in  air  parcels  confined  within  the 
Arctic  vortex.  They  showed  that  as  much  as  half  of  the 
vortex  air  could  be  processed  below  the  ice  frost  point 
by  December  20,  1994,  while  synoptic  temperatures 
never  fell  this  low.  This  rate  of  occurrence  of  mountain 
waves  with  temperatures  below  the  frost  point  suggests 
that  they  might  act  as  a  significant  source  of  localized 
ice  clouds  and  hence  solid  nitric  acid  hydrate  particles. 

Tabazadeh  et  al.  [1996]  have  examined  aircraft  lidar 
data  for  possible  influences  of  mesoscale  activity.  They 
concluded  that  only  about  half  of  the  PSCs  containing 
solid  particles  could  be  attributed  to  mountain  wave 
cooling  events.  However,  they  restricted  their  analysis 
to  a  single  day  and  to  a  relatively  small  region  of  the 
northern  hemisphere.  Here  we  consider  the  potential 
impact  of  mountain  waves  for  a  2  month  period  over  the 
entire  Arctic  vortex.  We  also  show  that  if  scenarios  4  or 
5  operate,  then  NAT  formation  may  occur  in  air  parcels 
even  several  weeks  after  an  initial  wave  encounter,  thus 
making  an  analysis  based  on  a  limited  number  of  back 
trajectories  very  difficult. 

2.  Method  of  Calculation 

We  calculate  the  frequency  with  which  mountain 
waves  are  encountered  by  air  parcels  confined  within 
the  Arctic  polar  vortex.  The  flow  of  air  parcels  was 
calculated  using  500  vortex-filling  trajectories  on  three 
different  potential  temperature  levels.  By  coupling  this 
information  with  several  realistic  NAH  particle  forma¬ 
tion  mechanisms  we  assess  the  potential  for  mountain 
waves  to  generate  solid  PSC  particles. 

2.1.  Mountain  Wave  Activity 

Calculation  of  mountain  wave  cooling  over  such  a 
large  geographical  area  and  for  a  whole  winter  first 
requires  a  reliable  parameterization  of  stratospheric 
mountain  wave  activity.  We  use  the  orographic  grav¬ 
ity  wave  parameterization  of  Bacmeister  et  al.  [1994]. 
This  is  based  on  a  global  database  of  dominant  topo¬ 
graphic  “ridge”  features,  which  contains  estimates  of 
ridge  height,  width,  and  orientation.  Synoptic  analy¬ 
ses  of  wind  and  temperature  at  standard  pressure  levels 
between  the  ground  and  10  mbar  are  then  used  to  deter¬ 
mine  air  flow  over,  these  ridge  features,  and  subsequent 
generation  of  mountain  waves.  Here  we  use  the  National 
Center  for  Environmental  Prediction  (NCEP)  analyses 
at  5°  x  2°  resolution.  A  two-dimensional  (2-D)  hydro¬ 
static  gravity  wave  model  is  used  in  conjunction  with 
the  NCEP  winds  and  temperatures  to  simulate  the  sub¬ 
sequent  amplification,  dissipation,  and/or  critical-level 
filtering  of  these  waves  as  they  propagate  away  from  the 
ridge  to  higher  regions  of  the  atmosphere. 

Advection  of  air  by  hydrostatic  gravity  waves  is  isen- 
tropic  to  a  good  approximation,  so  that  air  parcel  tern- 
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peratures  vary  in  response  to  wave-induced  vertical  ad- 
vection  according  to  the  dry  adiabatic  lapse  rate  [Eck- 
ermann  et  a/.,  1998].  At  any  instant,  this  yields  a  (re¬ 
versible)  temperature  perturbation  AT  from  the  back¬ 
ground  temperature.  In  the  present  study,  it  is  the  adi¬ 
abatic  cooling  (AT  <  0)  that  is  of  primary  significance, 
since  reversible  episodes  of  wave-induced  cooling  have 
major  effects  on  PSC  microphysics  and  chemistry.  We 
neglect  the  adiabatic  warming  phase  in  the  mountain 
waves  (AT  >  0)  as  a  loss  mechanism  of  nucleated  NAH 
particles.  This  warming  is  often  of  short  enough  dura¬ 
tion  in  mountain  waves  that  NAH  particles  do  not  com¬ 
pletely  evaporate  before  temperatures  fall  again  [Tsias 
et  a/,,  1997;  Carslaw  et  a/.,  1998a]. 

The  intrinsic  frequency  of  the  wave  u>  determines 
the  adiabatic  heating/cooling  rate  perturbation,  dT/dt , 
within  the  parcel.  However,  we  do  not  use  this  in¬ 
formation  to  simulate  scenario  2  since  the  high  heat¬ 
ing/cooling  rates  required  for  solid  particle  formation 
in  gravity  waves  [Tsias  et  a/.,  1997]  are  likely  to  occur 
under  non  hydrostatic  conditions  and  cannot  therefore 
be  described  by  this  simplified  model. 

The  maximum  wave-induced  adiabatic  cooling 
ATmax  (i.e.,  maximum  negative  departure  of  the  air 
parcel  temperature  from  that  expected  on  the  synoptic 
scale)  was  calculated  for  each  ridge  element  that  gener¬ 
ated  a  mountain  wave  which  reached  the  stratosphere 
on  a  given  day.  These  values  ATmax(xo,  z)  were  derived 
from  the  wave  amplitudes  in  the  mountain  wave  model, 
which  provides  a  vertical  profile  ( z  =  altitude)  centered 
directly  above  the  central  ridge  position  xo.  A  horizon¬ 
tal  “influence  function”  was  then  applied  to  simulate 
the  likely  horizontal  spread  of  mountain  wave  activity 
some  distance  away  from  the  ridge  center.  Here  we  used 

A T(x,  z)  =  ATmax(io,z)  exp  [-(a;  -  a;0)2/w2]  (1) 

to  simulate  cooling  downsteam  of  the  ridge  (x  >  xo), 
and 

A T(x,  z)  =  ATmax(i0,  z)  exp  [-10(x  -  x0)2/w2]  (2) 

for  cooling  values  upstream  of  the  ridge  (x  <  xo),  where 
x  is  the  downstream  distance  from  the  “ridge  line”  in 
the  orthogonal  or  “cross-ridge”  direction,  and  w  is  the 
cross-ridge  width.  This  Gaussian  influence  function  rep¬ 
resents  a  simplified  but  plausible  wave  amplitude  enve¬ 
lope  above  a  2D  ridge,  with  more  energy  found  down¬ 
stream  of  the  ridge  than  upstream,  as  suggested  by  case 
studies.  The  temperature  change  is  assumed  to  be  zero 
in  extension  of  the  ridge  axis. 

2.2.  Vortex-Filling  Trajectories 

The  technique  of  using  multiple  trajectories  that  fill 
a  region  of  the  atmosphere  (domain  filling)  provides  an 
ideal  means  of  assessing  PSC  formation  by  small  ridge 
elements.  We  have  used  500  isentropic  60-day  trajecto¬ 
ries  on  the  400,  475,  and  550  K  potential  temperature 


surfaces,  chosen  to  be  confined  within  the  Arctic  po¬ 
lar  vortex  between  December  1,  1994  and  January  31, 
1995.  Trajectories  were  calculated  from  NCEP  analy¬ 
ses  with  a  1  hour  time  step.  The  air  parcels  remained 
evenly  distributed  within  the  vortex  during  this  period, 
thus  providing  a  reliable  statistical  view  of  the  flow  over 
mountain  waves  on  a  given  potential  temperature  sur¬ 
face.  This  technique  correctly  describes  the  rate  of  air¬ 
flow  through  individual  mountain  waves  on  a  given  po¬ 
tential  temperature  level  as  well  as  correctly  describing 
the  lifetime  of  solid  particles  in  air  parcels  that  undergo 
synoptic  cooling  and  heating  cycles  as  they  circle  the 
pole. 

We  have  used  isentropic  trajectories  so  that  the  PSC 
formation  resulting  from  mountain  waves  can  be  con¬ 
veniently  identified  with  a  given  potential  temperature 
surface.  Clearly,  the  afr  flow  is  actually  diabatic  dur¬ 
ing  the  2  month  period  being  considered,  with  the  air 
sinking  perhaps  as  much  as  4  or  5  km  between  De¬ 
cember  1  and  January  31.  However,  the  assumption 
of  isentropic  flow  need  only  be  valid  for  the  short  time 
that  nitric  acid  hydrate  particles  exist  in  an  air  parcel, 
which,  due  to  periodic  warming  over  the  Alaskan  sec¬ 
tor  of  the  polar  vortex,  is  typically  less  than  2-3  days 
in  the  Arctic.  The  60-day  trajectories  can  therefore 
be  considered  equivalent  to  a  sequence  of  many  shorter 
quasi-isentropic  trajectories.  We  have  also  compared 
the  results  obtained  with  these  isentropic  trajectories 
with  similar  calculations  that  take  account  of  diabatic 
motion.  The  results  confirmed  the  validity  of  assuming 
short-term  quasi-isentropic  flow. 

The  mountain  wave  database  was  searched  for  coin¬ 
cidences  at  every  1-hour  trajectory  step  and  for  every 
trajectory.  A  coincidence  occurs  where  a  straight  line 
connecting  the  current  and  previous  trajectory  loca¬ 
tions  cuts  the  influence  region  of  one  or  more  waves. 
Where  two  or  more  waves  overlapped,  the  greater  of 
the  two  peak  temperature  changes  (ATmax)  was  taken. 
The  peak  temperature  change  was  linearly  interpolated 
between  pressure  levels  to  obtain  A Tmax  for  the  pres¬ 
sure  level  of  the  trajectory.  The  minimum  air  parcel 
temperature  between  two  trajectory  points  was  deter¬ 
mined  along  a  straight  line  cutting  the  mountain  wave 
influence  region.  The  formation  of  nitric  acid  trihydrate 
(NAT)  by  various  mechanisms  was  then  examined  at 
this  minimum  temperature.  We  have  made  similar  cal¬ 
culations  with  as  few  as  100  trajectories,  with  very  sim¬ 
ilar  results.  This  suggests  that  the  chosen  500  trajec¬ 
tories  provide  a  good  statistical  representation  of  solid 
particle  formation  in  mesoscale  cooling  events. 

2.3.  Formation  of  Solid  PSC  Particles 

The  formation  of  NAT  particles  in  the  mountain 
waves  was  investigated  based  on  scenarios  1,  4,  and  5. 
Ice  and  NAT  equilibrium  temperatures  were  calculated 
according  to  Hanson  and  Mauersberger  [1988]  assuming 
a  water  vapor  mixing  ratio  that  varied  with  altitude  as 
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follows:  400  K  (4.2  ppmv);  475  K  (4.6  ppmv);  550  K 
(5.4  ppmv).  In  the  absence  of  a  representative  HNO3 
altitude  profile  a  constant  10  ppbv  HNO3  on  all  levels 
was  used. 

According  to  scenario  1,  ice  was  assumed  to  nucle¬ 
ate  in  liquid  droplets  4  K  below  the  ice  frost  point 
(Tfrost  -  4  K)  allowing  NAT  to  form  heterogeneously. 
The  remaining  small  amount  of  H2SO4  and  H2O  in 
each  frozen  particle  was  assumed  to  remain  liquid,  lead¬ 
ing  to  the  release  of  NAT  particles  internally  mixed 
with  H2SO4/H2O  upon  ice  evaporation.  Such  parti¬ 
cles  would  be  composed  principally  of  NAT  and  would 
appear  in  lidar  observations  as  solid  particles  due  to 
their  depolarization  of  the  return  signal.  Upon  further 
warming  above  the  NAT  equilibrium  temperature,  pure 
liquid  droplets  are  released,  thus  resetting  the  particle 
phase  to  the  initial  conditions. 

Note  that  we  do  not  say  anything  about  the  num¬ 
ber  density  of  NAT  particles  in  a  given  air  parcel,  but 
only  whether  the  air  parcel  contains  such  particles  or 
not.  In  fact,  number  densities  of  NAT  particles  exit¬ 
ing  mountain  waves  are  known  to  be  quite  low  [Carslaw 
et  ai,  1998a].  We  also  do  not  say  anything  about  the 
size  of  the  individual  NAT  particles,  which  may  often 
be  far  from  thermodynamic  equilibrium  with  HNO3  in 
the  gas  phase  [Peter,  1997;  Tsias  et  ai,  1997;  Carslaw 
et  ai,  1998a].  The  effect  of  slow  growth  and  evapora¬ 
tion  kinetics  probably  acts  to  increase  the  geographical 
area  containing  solid  particles,  which  is  not  taken  into 
account  in  the  present  simulation. 

The  amount  of  supercooling  with  respect  to  ice  before 
crystalization  is  an  important  parameter  in  the  model 
since  this  represents  a  barrier  to  ice  and  NAT  formation 
in  the  waves,  thus  reducing  the  amount  of  NAT  that 
can  be  released  onto  the  synoptic  scale.  For  example, 
for  typical  conditions  on  the  50  mbar  level,  tempera¬ 
tures  must  fall  as  low  as  183  K  before  ice  formation, 
which  is  a  considerable  cooling  below  typical  synoptic 
temperatures  in  the  Arctic  winter  stratosphere.  Our  as¬ 
sumption  of  4  K  supercooling  is  in  agreement  with  the 
analysis  of  mountain  wave  PSCs  performed  by  Carslaw 
et  ai  [1998a]  but  larger  than  the  approximate  2  K  sug¬ 
gested  by  the  theoretical  work  of  Tabazadeh  et  ai  [1997] 
and  the  earlier  experiments  on  H2SO4/H2O  droplets 
[Bertram  et  ai,  1996].  More  recent  laboratory  experi¬ 
ments  examining  the  freezing  of  H2SO4/H2O  droplets 
suggest  that  ice  forms  approximately  3  K  below  the  ice 
frost  point  under  typical  stratospheric  conditions  [Koop 
et  ai,  1998],  a  somewhat  greater  supercooling  than  sug¬ 
gested  by  Bertram  et  ai  [1996].  We  therefore  consider 
our  assumed  4  K  supercooling  to  be  an  upper  limit  to 
the  supercooling  requirement.  In  the  sensitivity  calcu¬ 
lations  we  explore  the  effect  of  changing  this  quantity. 

According  to  scenarios  4  and  5,  SAT  particles  can 
serve  as  sites  for  NAT  nucleation  at  temperatures  a  few 
K  above  the  ice  frost  point.  To  simulate  these  mech¬ 
anisms,  we  assume  that  ice  formation  in  the  moun¬ 
tain  waves  leads  to  complete  crystallization  of  the  re¬ 


maining  liquid  HNO3-H2SO4-H2O  to  form  NAT  and 
SAT  (see  Figure  1).  Subsequent  warming  then  releases 
NAT/SAT  mixed  particles  above  Tfrost  and  pure  SAT 
particles  above  Tnat-  These  SAT  particles  can  then 
support  NAT  formation  according  to  either  scenario  4 
or  5.  Scenario  4  was  simulated  by  allowing  NAT  to 
form  at  a  supersaturation  of  10  (approximately  3  K 
below  Tnat)-  In  scenario  5,  NAT  forms  upon  cooling 
below  the  SAT  deliquescence  point,  which  is  typically 
3-4  K  below  Tnat  (the  SAT  deliquescence  temperature 
was  calculated  in  terms  of  the  HNO3  and  H20  partial 
pressures  using  the  thermodynamic  model  of  Carslaw 
et  ai  [1995]).  The  effect  of  these  mechanisms  is  to  am¬ 
plify  the  NAT  production  due  to  mountain  waves  alone; 
SAT  particles  that  are  released  from  mountain  waves 
and  suffer  a  further  cooling  cycle  can  serve  as  sites  for 
NAT  nucleation  at  several  degrees  higher  temperature 
than  would  otherwise  be  required  if  only  mechanism  1 
operated.  These  “additional”  NAT  particles  may  form 
either  on  the  synoptic  scale  when  temperatures  fall  suf¬ 
ficiently  low  or  in  mountain  waves  that  are  too  warm  to 
support  ice  nucleation. 

Air  parcels  containing  NAT  particles  were  advected 
downstream  by  the  domain-filling  trajectories  and  were 
assumed  to  persist  in  the  air  parcel  until  temperatures 
rose  above  Tnat>  at  which  point  NAT  was  assumed  to 
evaporate  immediately.  Some  mountain  waves  are  suf¬ 
ficiently  vigorous  that  ice  and  NAT  particles  may  form 
in  the  wave  although  synoptic  temperatures  are  higher 
than  the  NAT  equilibrium  temperature.  No  production 
of  persistent  synoptic  scale  NAT  occurs  in  such  cases, 
although  these  waves  may  still  contribute  to  chlorine 
activation  in  the  air  parcel  [Carslaw  et  ai,  1998b], 

3.  Mountain  PSC  Formation  and  NAT 
Production 

3.1.  Scenario  1 

Plate  1  shows  maps  of  NAT  occurrence  (black 
streaks)  for  December  1994  and  January  1995  together 
with  locations  where  temperatures  fell  below  Tfrost  —  4K 
(red  dots).  Figure  2a  shows  a  time  series  of  the  num¬ 
ber  of  model  air  parcels  containing  NAT  particles  as 
a  fraction  of  the  500  trajectories  on  the  475  K  level. 
The  distribution  of  ice-forming  mountain  waves  (Plate 
1),  and  hence  NAT  production,  is  strongly  biased  to¬ 
ward  the  European  and  Siberian  sector  of  the  northern 
hemisphere,  where  low  synoptic  temperatures  coincide 
with  strong  wave  activity.  Although  strong  mountain 
wave  cooling  occurs  over  Alaska,  synoptic  temperatures 
in  this  part  of  the  polar  vortex  (Aleutian  high)  are  per¬ 
sistently  warmer  than  over  the  European  sector,  which 
prevents  ice  formation  in  these  waves. 

Gravity  wave  amplitudes  often  increase  with  increas¬ 
ing  altitude  due  to  the  reduction  in  air  density,  hence 
the  number  of  very  cold  waves  (T  <  Tfrost  —  4  K)  and 
their  geographical  spread  increases  with  increasing  al- 
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titude  (potential  temperature).  On  the  400  K  level 
(~16  km  altitude),  most  of  these  events  occur  over  the 
Siberian  mountains  at  90°  —  100°E  with  a  few  on  the  east 
coast  of  Greenland.  On  the  475  K  level  (~19  km)  the 
mountains  of  East  Greenland,  Scandinavia,  the  Urals, 
and  Novaya  Zemlya  are  all  apparent,  while  on  the  550  K 
level  (~22  km)  strong  mountain  waves  are  apparent 
even  over  Baffin  Island  (60°  —  70°W)  and  farther  west. 

The  Scandinavian  mountains  act  as  a  significant 
source  of  NAT  particles  that  extend  over  Scandinavia, 
particularly  on  the  475  K  level.  The  east  coast  of  Green¬ 
land  also  generates  NAT  particles  that  survive  the  pas¬ 
sage  over  the  North  Atlantic,  where  temperatures  are 
often  below  Tnat  over  a  wide  area.  PSCs  containing 
solid  particles  observed  over  Scandinavia  may  therefore 
often  originate  several  thousand  kilometers  upstream. 
NAT  particles  generated  over  the  Siberian  mountains 
(90°  E)  generally  persist  only  a  few  hours  before  evap¬ 
orating  in  the  warmer  temperatures  that  are  persistent 
over  the  Alaskan  sector  of  the  polar  vortex.  Although 
the  wave  amplitudes  and  number  of  very  cold  waves 
are  greatest  on  the  550  K  level,  the  number  of  model 
air  parcels  containing  NAT  is  very  small,  reflecting  the 
generally  higher  synoptic  temperatures  on  that  level  rel¬ 
ative  to  the  NAT  equilibrium  temperature.  The  fewer 
number  of  air  parcels  containing  NAT  on  the  400  and 
550  K  levels  is  also  due  to  the  location  of  the  polar  vor¬ 
tex  relative  to  significant  orography.  On  the  475  K  level 
the  area  where  temperatures  below  Tnat  occur  covered 
most  of  the  region  between  60°W  —  150°E  and  north 
of  65° N.  In  contrast,  the  potential  NAT  area  on  the 
550  K  level  was  restricted  to  ~  20°  W  (the  East  Green¬ 
land  coast),  110°E  and  north  of  ~  68°N.  Hence,  on 
the  550  K  level,  NAT  particle  persistence  downwind  of 
the  East  Greenland  coastal  mountains  and  the  northern 
Scandinavian  mountains  was  severely  restricted. 

The  number  of  air  parcels  containing  NAT  (Fig¬ 
ure  2a)  fluctuates  in  response  to  the  extent  of  low  tem¬ 
peratures  in  the  vortex.  During  cold  periods,  both  the 
production  rate  of  NAT  in  the  mountain  waves  and  the 
lifetime  of  the  generated  NAT  particles  are  greater.  The 


Plate  1.  Geographical  distribution  of  mountain  wave- 
induced  NAT  PSCs  during  December  1994  and  January 
1995  on  the  levels  (a)  400  K,  (b)  475  K,  and  (c)  550  K 
according  to  scenario  1.  Air  parcels  calculated  to  con¬ 
tain  NAT  particles  are  marked  with  a  black  dot  at  each 
trajectory  step  (1  hour  intervals).  NAT  particles  were 
assumed  to  be  produced  in  air  parcels  where  the  tem¬ 
perature  passed  below  Tfrost  -  4  K  (red  dots)  leading  to 
ice  formation.  Note  that  the  red  dots  are  placed  at  the 
nearest  trajectory  point  location  and  therefore  may  be 
displaced  away  from  the  ridge  element  that  induced  the 
adiabatic  cooling  event.  Synoptic  temperatures  during 
the  2  month  period  were  such  that  no  NAT  formation 
was  assumed  to  occur  other  than  in  mountain  waves. 
The  calculation  refers  to  500  trajectories. 
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Figure  2.  (a)  Fraction  of  the  500  model  air  parcels  con¬ 
taining  NAT  particles,  plotted  at  1  hour  intervals.  Thin 
solid  line,  assuming  ice  forms  at  Tfrost  —  4  K  (scenario  1, 
corresponding  to  the  map  in  Plate  1);  dotted  line,  as¬ 
suming  NAT  forms  upon  SAT  deliquescence  (scenario 
5);  thick  solid  line,  a  sensitivity  test  using  scenario  1 
in  which  all  wave  amplitudes  were  reduced  by  a  factor 
of  2;  dashed  line,  a  sensitivity  test  using  scenario  1  but 
with  all  wave  amplitudes  reduced  by  a  factor  of  2  and 
synoptic  temperatures  reduced  by  2  K  overall,  (b)  Frac¬ 
tion  of  air  parcels  containing  SAT  particles  according  to 
scenario  5. 

fraction  of  model  air  parcels  containing  NAT  (induced 
through  the  action  of  mountain  waves)  is  as  much  as  5% 
on  the  475  K  level,  equivalent  to  a  total  area  of  6  x  105 
km2  if  these  air  parcels  are  considered  representative 
of  5%  of  the  vortex  air.  Much  lower  coverages  are  pre¬ 
dicted  for  the  400  and  550  K  levels.  If  NAT  were  allowed 
to  form  whenever  temperatures  fell  below  the  NAT  equi¬ 
librium  temperature,  then  the  maximum  areal  coverage 
of  NAT  in  the  vortex  would  be  about  20%.  This  implies 
that  NAT  is  present  in  approximately  25%  of  the  vortex 
area  where  NAT  is  stable.  Locally,  the  areal  coverage 
of  NAT  can  be  much  higher.  For  example,  downwind 
of  the  Scandinavian  mountains  more  than  50%  of  the 
model  air  parcels  are  calculated  to  contain  NAT  parti¬ 
cles  on  January  11,  1995  (day  42  in  Figure  2). 

3.2.  Scenarios  4  and  5 

The  effect  of  NAT  production  through  SAT  deliques¬ 
cence  (scenario  5)  leads  to  a  considerable  increase  in 


the  area  of  the  vortex  containing  NAT  (dotted  line  in 
Figure  2a);  about  a  factor  of  2  increase  compared  with 
scenario  1  in  mid- January,  and  sometimes  as  much  as  a 
factor  of  6  increase.  NAT  is  predicted  to  fill  up  to  14% 
of  the  vortex  on  the  475  K  level,  or  about  2  x  107  km2. 

It  is  interesting  to  note  that  the  amount  of  air  con¬ 
taining  SAT  particles  increases  steadily  during  Decem¬ 
ber  and  early  January  (Figure  2b).  Temperatures  dur¬ 
ing  this  period  generally  remain  below  the  SAT  melt¬ 
ing  temperature  of  210-215  K,  thus  SAT  generated  in 
mountain  waves  accumulates  in  the  vortex  (note  that 
NAT  formation  upon  SAT  deliquescence  does  not  lead 
to  the  loss  of  SAT  particles  since  these  are  assumed  to 
be  recovered  once  the  NAT  evaporates  again).  By  mid- 
January  as  much  as  50%  of  the  vortex  is  calculated  to 
contain  some  SAT.  The  amount  of  SAT  falls  only  in  late 
January  as  temperatures  increase  and  the  loss  rate  of 
SAT  by  melting  exceeds  the  production  rate  in  moun¬ 
tain  waves.  Note  that  the  assumption  of  isentropic  air 
parcel  behavior  is  not  strictly  valid  for  the  much  longer 
times  (up  to  several  weeks)  that  SAT  particles  exist  in 
a  given  air  parcel.  To  test  the  results  in  Figure  2b,  we 
have  performed  similar  calculations  based  on  diabatic 
trajectories  cutting  the  475  K  level  (as  used  by  Carslaw 
et  al.  [1998b])  and  obtained  very  similar  results  to  those 
shown  here. 

Particularly  important  in  terms  of  being  able  to  at¬ 
tribute  NAT  formation  in  the  stratosphere  to  a  partic¬ 
ular  mechanism  is  that  in  many  of  the  air  parcels  NAT 
formation  upon  SAT  deliquescence  occurred  several 
weeks  after  the  initial  mountain  wave  event.  Attribut¬ 
ing  a  given  solid-particle  PSC  observation  to  mountain 
waves  might  therefore  prove  very  difficult  if  this  SAT- 
deliquescence  mechanism  operates.  Such  NAT  PSCs 
may  also  be  more  uniform  in  structure  than  those 
derived  directly  from  mountain  waves  due  to  the  at¬ 
mospheric  mixing  that  occurs  during  the  few  weeks. 

An  alternative  to  the  SAT  deliquescence  mechanism 
is  the  preactivation  mechanism  of  Zhang  et  al.  [1996] 
(scenario  4).  In  their  mechanism,  NAT  can  form  on 
SAT  at  a  supersaturation  with  respect  to  NAT,  Snat,  of 
~10.  This  corresponds  to  a  temperature  of  about  2-3  K 
below  the  NAT  equilibrium  temperature,  Ts at?  under 
normal  stratospheric  conditions.  Formation  of  NAT  by 
this  process  is  therefore  likely  to  occur  more  often  than 
by  the  SAT  deliquescence  mechanism,  which  requires 
that  temperatures  fall  about  4  K  below  Tnat-  Whether 
NAT  formation  due  to  SAT  preactivation  or  SAT  del¬ 
iquescence  operates  in  the  stratosphere  is  not  certain. 
A  sensitivity  test  in  which  NAT  formed  via  the  pre¬ 
activation  mechanism  led  to  calculated  NAT  amounts 
generally  about  30%  higher  than  predicted  by  assum¬ 
ing  NAT  formation  upon  SAT  deliquescence  (20%  NAT 
coverage  of  the  vortex  compared  with  15%). 

3.3.  Sensitivity  Tests 

There  are  several  factors  that  may  change  the  produc¬ 
tion  of  solid  particles  initiated  by  mountain  waves.  We 
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show  how  these  might  affect  the  foregoing  calculations 
for  the  475  K  level  also  in  Figure  2a. 

3.3.1.  Reduction  in  wave  amplitudes.  We 
have  considered  the  effect  of  reducing  the  magnitude 
of  adiabatic  cooling  in  all  waves  (thick  solid  line  in  Fig¬ 
ure  2a).  Bacmeister  et  al.  [1994]  and  later  Carslaw 
et  al.  [1998a]  have  compared  the  predictions  of  the 
Bacmeister  et  al.  model  with  observations  of  moun¬ 
tain  wave  cooling.  They  found  reasonable  agreement 
between  predicted  and  observed  adiabatic  temperature 
excursions  for  large-amplitude  waves  where  AT  was  10- 
20  K.  However,  there  is  some  indication  that  the  model 
somewhat  overestimates  the  wave  amplitude  for  less 
vigorous  waves.  A  comparison  with  waves  observed  by 
the  ER-2  over  the  southern  tip  of  Greenland  suggests 
that  the  Bacmeister  et  al.  model  overestimated  the 
wave  amplitude  by  about  50%  (AT  ~  12  K  compared 
with  an  observed  8  K).  As  a  sensitivity  test,  we  have  re¬ 
duced  all  wave  amplitudes  by  a  factor  of  2.  This  leads 
to  a  considerable  reduction  in  the  amount  of  mountain 
wave-initiated  NAT.  NAT  is  predicted  to  occur  on  only 
a  few  days  in  January  and  to  be  present  in  less  than 
~1.5%  of  the  vortex,  or  1.5  x  105  km2.  The  maximum 
spatial  extent  of  SAT  particles  in  the  vortex  reaches 
about  6%  of  the  vortex  area,  hence  even  allowing  for 
NAT  production  through  SAT  deliquescence  does  not 
significantly  increase  the  amount  of  NAT. 

3.3.2.  Bias  in  synoptic  temperatures.  The 
dashed  line  in  Figure  2a  shows  the  results  of  a  simu¬ 
lation  with  all  mountain  wave  amplitudes  reduced  by  a 
factor  of  2,  as  above,  but  with  the  synoptic  tempera¬ 
tures  decreased  by  2  K  overall.  This  leads  to  NAT  cov¬ 
erage  in  the  vortex  as  much  as  a  factor  of  3-4  greater 
than  predicted  by  reducing  the  wave  amplitudes  alone 
and  reveals  the  enormous  sensitivity  of  NAT  produc¬ 
tion  to  temperature,  which  arises  due  to  the  combined 
increase  in  the  number  of  ice  PSCs  occurring  in  the 
waves  as  well  as  the  extended  persistence  of  NAT  on 
the  synoptic  scale. 

3.3.3.  Ice  supercooling  requirement.  A  fur¬ 
ther  test  was  made  by  reducing  the  amount  of  super¬ 
cooling  required  for  ice  formation  in  the  waves  from  4  K 
to  2  K  (results  not  shown).  Changing  the  supercooling 
requirement  for  ice  formation  affects  only  the  efficiency 
of  NAT  production  but  not  the  persistence  of  NAT  on 
the  synoptic  scale.  The  fraction  of  model  air  parcels 
containing  NAT  was  therefore  between  that  predicted 
for  the  sensitivity  tests  described  in  3.3.1  and  3.3.2.  In¬ 
creasing  the  amount  of  water  vapor  on  the  475  K  level 
to  5  ppmv  changes  the  ice  frost  point  and  NAT  equi¬ 
librium  temperatures  by  less  than  0.5  K  and  therefore 
does  not  significantly  affect  these  results. 

4.  Discussion 

The  model  simulations  presented  above  reveal  that 
NAT  production  due  to  mountain  waves  depends  very 
sensitively  on  a  number  of  parameters.  The  most  sen¬ 
sitive,  in  order  of  importance,  are  (1)  the  magnitude  of 


adiabatic  cooling  in  each  wave,  (2)  the  synoptic  temper¬ 
ature,  and  (3)  the  supercooling  required  for  initial  ice 
formation.  The  model  simulations  suggest  that  some 
mountain  wave-induced  solid  particles  are  likely  to  be 
present  in  the  polar  vortex  when  temperatures  are  low 
enough  to  support  nitric  acid  hydrates,  but  the  current 
uncertainties  in  model  parameters  do  not  allow  the  cal¬ 
culated  areal  coverages  to  be  well  constrained  (cover¬ 
ages  of  NAT  are  calculated  to  lie  between  1  and  14%  of 
the  vortex  area  on  the  475  K  level,  or,  equivalently, 
between  about  4  and  60%  of  the  vortex  area  where 
temperatures  are  below  the  NAT  equilibrium  temper¬ 
ature).  Directly  downwind  of  the  Scandinavian  moun¬ 
tains,  more  than  50%  of  the  air  on  a  single  level  may 
contain  mountain  wave-induced  NAT  under  favorable 
conditions.  Our  calculations  did  not  take  account  of 
mixing  of  NAT- containing  air  parcels  into  the  surround¬ 
ing  air.  This  would  serve  to  increase  the  fraction  of  the 
polar  vortex  containing  some  NAT  particles. 

The  magnitude  of  the  adiabatic  cooling  is  a  par¬ 
ticularly  uncertain  parameter.  Our  test  calculations 
assumed  that  the  adiabatic  cooling  predicted  by  the 
Bacmeister  et  al.  model  might  be  too  large  by  up  to 
a  factor  of  2.  This  is  equivalent  to  a  reduction  in  the 
cooling  in  large  waves  from  AT  ~  20K  to  AT  ~  10K, 
a  far  greater  uncertainty  than  in  the  supercooling  re¬ 
quired  for  initial  ice  formation  (AT  ~  4K).  The  effect 
of  this  reduction  in  model  wave  amplitudes  was  to  re¬ 
duce  the  areal  coverage  of  NAT  by  more  than  a  fac¬ 
tor  of  3,  or  even  to  switch  off  NAT  production  com¬ 
pletely.  Adiabatic  cooling  in  mountain  waves  by  as 
much  as  AT  ~  12  —  19K  has  been  inferred  from  ob¬ 
servations  of  ice  clouds  over  Norway  [Carslaw  et  a/., 
1998a,b],  suggesting  that  a  uniform  reduction  in  pre¬ 
dicted  wave  amplitudes  may  be  unrealistic.  Further 
work  is  clearly  required  to  validate  the  mountain  wave 
forecast  model  that  we  have  used.  However,  it  is  worth 
pointing  out  that  a  realistic  reduction  in  model  wave 
amplitudes,  which  would  result  by  introducing  disper¬ 
sion  of  wave  energy  in  the  horizontal,  will  lead  to  a  much 
wider  influence  region  around  each  ridge  element  in  the 
model.  This  in  turn  will  result  in  more  vortex  air  being 
processed  through  mountain  waves,  which  might  partly 
compensate  for  the  reduced  peak  cooling. 

There  is  clear  evidence  for  the  formation  of  ice  in 
mountain  waves  and  the  subsequent  release  of  solid 
particles  downwind  of  the  mountains  [Carslaw  et  a/., 
1998a].  However,  the  crystallization  of  the  remaining 
H2SO4  to  form  sulfuric  acid  hydrates  such  as  SAT  is  less 
certain.  While  SAT  formation  in  mountain  waves  is  not 
a  requirement  for  solid  particle  formation,  we  showed 
that  subsequent  deliquescence  of  SAT  upon  cooling  (fol¬ 
lowed  by  NAT  formation  [Koop  and  Carslaw ,  1996])  or 
nucleation  of  NAT  on  preactivated  SAT  particles  [Zhang 
et  a/.,  1996]  could  act  as  an  efficient  additional  NAT 
source. 

We  have  also  highlighted  that  the  fraction  of  the  polar 
vortex  containing  mountain  wave-generated  NAT  par- 
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tides  depends  on  the  location  of  low  temperatures  rel¬ 
ative  to  the  Greenlandic  and  Scandinavian  mountains. 
In  particular,  the  east  coast  of  Greenland  can  serve  as 
a  source  of  solid  PSC  particles  that  flow  over  Europe 
provided  that  temperatures  are  below  Tnat  over  the 
North  Atlantic.  Interannual  changes  in  the  location  of 
the  vortex  and  associated  low  temperatures  could  there¬ 
fore  strongly  influence  the  amount  of  mountain  wave- 
generated  solid-containing  PSCs. 

A  related  factor  is  that  the  largest  amplitude  gravity 
waves  are  sustained  in  the  stratosphere  under  condi¬ 
tions  where  the  stratospheric  winds  are  strongest.  The 
location  of  the  polar-night  jet  (or  polar  vortex  edge) 
relative  to  the  region  of  lowest  stratospheric  synoptic 
tempertures  is  therefore  also  important.  Particularly 
important  for  the  effectiveness  of  mountain  waves  in 
generating  NAT  is  that  the  coldest  region  of  the  po¬ 
lar  vortex  can  sometimes  lie  close  to  or  on  the  vortex 
edge.  Under  these  conditions,  the  potential  for  sustain¬ 
ing  large  wave  amplitudes  (highest  stratospheric  winds) 
combines  with  the  maximum  likelihood  that  NAT  will 
presist  on  the  synoptic  scale  (lowest  synoptic  tempera¬ 
tures).  Considerable  NAT  production  due  to  mountain 
waves  might  therefore  be  expected  during  these  periods. 

Direct  comparison  between  observations  of  solid- 
containing  PSCs  and  the  predicted  location  of  NAT 
clouds  originating  from  mountain  waves  should  be  made 
with  caution.  The  mountain  wave  climatology  used  here 
should  predict  the  general  location  and  magnitude  of 
wave  activity  on  a  statistical  basis.  However,  individ¬ 
ual  back  trajectories  constructed  for  comparison  with 
PSC  observations  cannot  be  expected  to  yield  reliable 
statistical  information  on  NAT  production.  We  have 
also  shown  that  NAT  formation  as  a  result  of  SAT  del¬ 
iquescence  close  to  the  frost  point  [Koop  and  Carslaw, 
1996;  Iraci  et  al.,  1998]  could  disguise  any  direct  link  be¬ 
tween  observed  solid  particles  and  the  initiating  moun¬ 
tain  wave. 

Larsen  et  ai  [1997]  have  performed  the  most  com¬ 
plete  statistical  analysis  of  the  synoptic-  conditions  un¬ 
der  which  PSCs  are  likely  to  be  observed  in  the  solid  or 
liquid  phase.  Their  investigation,  based  on  30  balloon- 
sonde  flights  over  4  years  from  seven  ground  stations  in 
the  Arctic,  reveals  that  the  likelihood  of  solid  and  liq¬ 
uid  PSCs  depends  very  much  on  the  synoptic  conditions 
prior  to  the  observation.  In  particular,  the  occurrence 
of  PSCs  containing  solid  particles  was  shown  to  be  more 
likely  in  air  parcels  that  had  spent  a  considerable  time 
below  the  NAT  equilibrium  temperature,  while  liquid 
aerosols  occur  most  frequently  in  air  parcels  that  have 
cooled  recently.  Our  proposed  mechanism,  in  terms  of 
mountain  wave  NAT  activation,  is  qualitatively  consis¬ 
tent  with  these  statistics;  the  more  time  air  spends  at 
low  temperatures  the  greater  is  the  probability  that  a 
mountain  wave  will  force  the  temperature  low  enough 
to  induce  ice  formation  and  hence  solidification  of  the 
HNO3-H2SO4-H2O  solutions.  The  likelihood  of  NAT 
activation  upon  SAT  deliquescence  is  also  greater  at 


lower  temperatures  since  the  deliquescence  temperature 
is  about  3-5  K  below  the  NAT  equilibrium  temperature. 

Future  work  should  extend  the  study  presented  here 
to  include  a  quantitative  comparison  with  PSC  obser¬ 
vations.  Such  a  study  should  aim  to  establish,  as  far  as 
possible,  the  number  of  solid  particle-containing  PSC 
observations  that  can  be  explained  by  mountain  wave 
processes.  One  approach  would  be  to  compile  a  lidar 
climatology,  with  PSCs  being  classified  according  to  a 
threshold  value  of  depolarization. 
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